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Abstract. This investigation studies the decidability problem of plane edge 
coloring with three symbols. In the edge coloring (or Wang tiles) of a plane, 
unit squares with colored edges that have one of p colors are arranged side by 
side such that the touching edges of the adjacent tiles have the same colors. 
Given a basic set B of Wang tiles, the decision problem is to find an algorithm 
to determine whether or not S(-B) ^ 0, where is the set of all global 

patterns on Ip that can be constructed from the Wang tiles in B. 

When p > 5, the problem is known to be undecidable. When p = 2, the 
problem is decidable. This study proves that when p = 3, the problem is 
also decidable. ViB) is the set of all periodic patterns on Z-^ that can be 
generated by the tiles in B. If P(-B) 0, then B has a subset B' of minimal 
cycle generators such that •p(B') ^ and ViB") = for B" g B' . This 
study demonstrates that the set C(3) of all minimal cycle generators contains 
787, 605 members that can be classified into 2, 906 equivalence classes. J\f{3) is 
the set of all maximal non-cycle generators : if _B £ A''(3), then V{B) = and 
V{B) for B ^ B. The problem is shown to be decidable by proving that 
B G A/'(3) impUes = 0. Consequently, E(B) if and only if 7^ 0. 



1. Introduction 

Decidability problems have been studied for many years. See, for example, the 
recent review article of Goodman-Strauss ^ . One of the most active areas of re- 
search into the decidability problem is that of the plane tiling [7] . This study focuses 
on the decidability problem concerning plane edge coloring with three symbols. 

The coloring of using unit squares has a long history [7]. In 1961, Wang 
jT4] started to study the square tiling of a plane to prove theorems by pattern 
recognition. Unit squares with colored edges are arranged side by side so that the 
touching edges of the adjacent tiles have the same color; the tiles cannot be rotated 
or reflected. Today, such tiles are called Wang tiles or Wang dominos [H [7] . 

The 2x2 unit square is denoted by Z2x2- The set of p colors is {0, 1, ■ • • ,p—l}. 
Therefore, the total set of Wang tiles is denoted by i;2x2(p) = {0, 1, • • • ,p- l}^^x^ 
A set B of Wang tiles is called a basic set (of Wang tiles). Let T,{B) and V{B) be 
the sets of all global patterns and periodic patterns on Z^, respectively, that can 
be constructed from the Wang tiles in B. 

The decision problem concerning tiling with of Wang tiles concerns the existence 
of an algorithm that can determine whether or not 

(1.1) S(B)^0 

for any finite set B of Wang tiles. 
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Clearly, ViB) C T.{B), meaning that iiP{B) ^ 0, then ^ 0. In [M], Wang 

conjectured that any set of tiles that can tile a plane can tile the plane periodically: 

(1.2) if ^ 0, then P{B) ^ 0. 

If (|1.2p holds, then the decision problem that is spcified by (jl.ip is reduced to 
the much easier problem of determining whether or not 

(1.3) ViB) ^9. 

However, in 1966, Berger [4] proved that Wang's conjecture was wrong and the 
decision problem concerning Wang's tiling is undecidable. He presented a set B of 
20426 Wang tiles that could only tile the plane aperiodically: 

(1.4) S(B) ^ and P{B) = 0. 

Later, he reduced the number of tiles to 104. Thereafter, smaller basic sets were 
found by Knuth, Lauchli, Robinson, Penrose, Ammann, Culik and Kari. Currently, 
the smallest number of tiles that can tile a plane aperiodically is 13, with five colors: 
(fO]) holds and then ([^2]) fails for p = 5 [5]. 

Recently, Hu and Lin f8^ showed that Wang's conjecture (|1.2|) holds ifp — 2: any 
set of Wang tiles with two colors that can tile a plane can tile the plane periodically. 

In that study, they showed that statement p.2p can be approached by studying 
how periodic patterns can be generated from a given basic set. First, B is called a 
cycle generator if ViB) ^ 0; otherwise, B is called a non-cycle generator. Moreover, 
B C S2x2(p) is called a minimal cycle generator (MCG) if i? is a cycle generator 
and ViB') = whenever B' C B\ B C E2x2(p) is called a maximal non-cycle 
generator (MNCG) if i? is a non-cycle generator and V{B") ^ for any B" ^ B. 

Given p > 2, denote the set of all minimal cycle generators by C{p) and the set 
of maximal non-cycle generators by Afip) . Clearly, 

(1.5) C(p)nA/'(p) = 0. 
Statement (fOj) follows for p > 2 if 

(1.6) E(B) ^ for any B e M{p) 

can be shown. Indeed, in ^8^, it is shown that C(2) has 38 members ;A/'(2) has nine 
members, and (|1.6p holds for p = 2. This paper studies the case of p = 3. Now, 
C(3) and A/'(3) have close to a million members and cannot be handled manually. 
After the symmetry group -D4 of ^2x2 and the permutation group Sp of colors of 
horizontal and vertical edges, respectively, are applied, C(3) still contains thousands 
of equivalent classes. Hence, computer programs are utilized to determine C(3) and 
Ar(3) and finally (|1.6I) is shown to hold for p = 3. Therefore, the problem (II. ip is 
decidable for p = 3. 

For p = 4, C(4) is enormous. Therefore, the arguments and the computer pro- 
gram need to be much efficient to handle this situation. 

Corner coloring with p — i can be treated similarly. The result will be presented 
elsewhere. 

The rest of paper is arranged as follows. Section 2 introduces the ordering 
matrix of all 81 local patterns and classifies them into three groups. The recurrence 
formula for patterns on Smxn are derived. It is important in proving (|1.6|) - that 
the maximum non-cycle generators cannot generate global patterns. Section 3 will 
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introduce the procedure for determining the sets C(3) and A/'(3). The main result 
is proven using a computer. 

2. Preliminary 

This section introduces all necessary elements for proving ()1.6p . First, let I]„ixn(-B) 
be the set of all local patterns on Zmxn that can be generated by B. Clearly, 

(2.1) if I]™xn(S) = for some m, n > 2, then = 0. 

2.1. Symmetries. The symmetry of the unit square ^2x2 is introduced. The 
symmetry group of the rectangle 7^2x2 is D^, which is the dihedral group of order 
eight. The group D4 is generated by rotation p through ^ and reflection m about 
the y-axis. Denote the elements of D4 by D4 = {I , , ,m,mp,mp^ ,mp^^ . 




Figure 2.1 

Therefore, given a basic set B C S2x2(p) and any element r G -D4, another basic 
set {B)t can be obtained by transforming the local patterns in B by r. 

Additionally, consider the permutation group S'p on {0, 1, • • • ,p — 1}. If 77 £ S'p 
and ry(0) = ip, ?7(1) = ii, ■ • • , ri(p — 1) = ip-i, we write 

/ 1 ... ,-iy 

For rj £ Sp and B 6 T,2x2{p), another basic set {B)^ can be obtained. 

In edge coloring, the permutations of colors in the horizontal and vertical direc- 
tions are mutually independent. Denote the permutations of colors in the horizontal 
and vertical edges by rjh G Sp and rjv G Sp, respectively. Now, for any B C ^2x2{p), 
define the equivalence class [B] of B by 

(2.2) [B] = {B' C S2x2(p) : B' = (((i?),) Vh)rh^''' ^ D4 and r]h,riy e Sp}. 

In [5], whether or not S(i?) = and P{B) = is shown to be independent of 
the choice of elements in [B]. Indeed, for any B' E [B], 

E(S') 7^ (or P(S') ^ 0) if and only if I](B) ^ (or 7^(5) ^ 0). 

Moreover, for B' e [B], B' is an MCG (MNCG) if and only if B is an MCG 
(MNCG). Therefore, groups D4 and S3 can be used efficiently to reduce the number 
of cases B C 1^2x2(3) that must be considered, greatly reducing the computation 
time. 



4 



HUNG-HSUN CHEN, WEN-GUEI HU, DE-JAN LAI, AND SONG-SUN LIN* 



2.2. Ordering Matrix. Now, the case p = 3 is considered. The vertical ordering 
matrix Y2x2 = [yi,j]gx9 local patterns in S2x2(p) is given by 



(2.3) 



■ 2x2 



M m 



m 



m 
m 



m 



m 



m m 
m m 



m m 
m m 



m m 



Q\l I/O' 



(2.4) 



Y2;l 


Y2;2 


Y2;3 


Y2;4 


Y2;5 


Y2;6 


Y2;7 


Y2;8 


Y2;9 



3x3 



The recurrence relation of Y^+i is easily obtained as follows. Denote by 

9 

(2.5) Y2 = ^ Y2;i 



and 

(2.6) 

where 



Y2;i — [ V2-i:,p,q ] 



3x3 ■ 




and i = 1 + ai • 3^ + 02 • 3", G {0, 1, 2}. For m > 2, denote by 



(2.7) 

and 
(2.8) 



Ym+l — Y,„+l;z 

1=1 
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where ym+i;i:p,q is the set of all patterns of the form 



\<71 / 


\ «2 / 














/ Pl\ 


/ P2\ 







where a ,/3, pfc and G {0, 1, 2}, 1 < fc < to, such that 

i = 1 + a ■ 31 + ^ • 3"^ 



(2.9) 



P = 1 + EPfc3 
fe=i 

m 

<Z = 1 + E 9/03' 



and • G {0, 1, 2}. Therefore, when i = 1, 4, 7, 



Y. 



3 

IZ y2;3 + i-l;l,lY„ 


i;3j-2 


3 


-1;1,2Y„ 


t;33-2 


3 

H a2:3 + j 


-1;1,3Y„ 


,;3j-2 


j=l 
3 




3 = 1 
3 






3 = 1 
3 






51 y2;3 + i-l;2.lY„ 


i;3j-2 


T, y2■.j+^ 


-1;2,2 Y„ 


t:33-2 


E a2;3 + ! 


-l;2,3Y'n 


i;33-2 


3 = 1 




3 = 1 






3 = 1 






3 




3 






3 






5D y2!3 + i-l;3,lY, 


i;33-2 




-1;3,2 Yn 


i;33-2 


E J/2;3 + . 


-1;3,3Y„ 


i;33-2 


J = l 




i=i 






J = l 







3^x3" 



when i = 2, 5, i 



3 3 
Yl y2;3 + i-2;l.lYm;33_i V 2 ;3 +i - 2: 1 , 2 Y^ jSj _ i 

3 = 1 3 = 1 

3 



Yl y2;3 + i-2;l,3Ym;3j_i 
3 = 1 



Y y2;3 + i-2;2.lYm;33-l J/2 ;3 + i - 2 ; 2 , 2 Y™ ; 3 j _ 1 J2 a2;3 + i - 2 ; 2 , 3 Y,„ ;3 j - 



3 = 1 3 = 1 

3 3 

Y y2i3 + i-2;3,lY„;33_i X] 2/2 ;3 + i - 2 : 3 , 2 Y„ ; 33 _ 1 

3 = 1 3 = 1 



Y »/2;3 + i-2;3,3Y„;3j_i 
3 = 1 



when i = 3, 6, 9, 



■ m+l;' 



3 

E y2;3 + i-3;l.lY„ 


i;33 


3 

Y y2;j + i~3;l,2Y„ 


i;33 


3 

Y i/2;j + i-3;l,3Y„ 


^;33 


3 = 1 




3 = 1 




3 = 1 




3 




3 




3 




Y y2;3 + i-3;2,lY„ 


i;33 


Y !/2;3+i~3;2,2 Y„ 


i;33 


Y !y2;3+i-3;2,3 Y„ 


^;33 


3 = 1 




3 = 1 




3 = 1 




3 




3 




3 




Y y2;3 + i-3;3,l Yt, 


i;33 


Y y2;3+i-3;3,2 Y„ 


i;33 


Y y2;3+i-3;3,3 Y„ 


x;33 


L 3 = 1 




3 = 1 




3 = 1 





Given B C S2x2(3), the associated vertical transition matrix V2x2(S) is defined 
by V2x2(-B) = [vi^j], where = 1 if and only if yij G B. 

The recurrence formula for a higher-order vertical transition matrix can be ob- 
tained as follows. Denote by 



9 

V2 - Y.^2■,^, 

i=l 



with 



For m > 2, denote by 



9 J3X3- 



9 

= y~^"V"m+l;i- 

i=l 
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Now, for i = 1,4, 7, 



m+l:', 



J2 f2;j + i-l 
j = l 
3 

E ■"2:j+i-l 
j = l 



J=l 

3 



.,2Vm;3j-2 E 't'2;i + i - 1 ; 1 , 3 ;3 j - 2 



2,lV„^.3j_2 E ^2;j + i-l;2,2 V^.3j_2 
J = l 



E ^2;j + i-l;3,l VT^;3j_2 E ^2 ; j + i - 1 ;3 , 2 V^^i ; 3j _ 2 



3 = 1 
3 

'"2;3 + i-l;2,3 Vm;3j_2 

3 = 1 
3 

Yl ''2;3 + s-l;3,3V„;3j_2 
3 = 1 



3'"x3" 



for i = 2,5, 8, 



for z = 3, 6, 9, 



E ''2;j-l-. 


-2,1. lV„ 


i;33-l 


E ''2;3 + . 


-2;1,2 


^;33-i 


E ''2;3 + . 


-2;1,3 


i;33-i 


j = l 






3 = 1 






3 = 1 






3 






3 






3 






E V2-j + , 


-2;2,1 V„ 


i;33-l 


E l'2;j + . 


-2;2,2 


^;33-i 


E ^2-.j + z 


-2;2,3 Vj 


^;33-i 


j = l 






3 = 1 






3 = 1 






3 






3 






3 






E "2;j-|- = 


-2;3,1 V„ 


i:33-l 


E V2;j + , 


-2;3,2 Vt, 


a;33-l 


E ^2;j + , 


-2;3,3 


a;33-l 


3 = 1 






3 = 1 






3 = 1 







3'"x3" 



r 3 3 3 

E "2;3-|-i-3!l,lV,„:3j E '^2:3 + , - 3; 1 ,2 V,„:3j E ^'2;3 + i - 3; 1 ,3 V„;3j 

3 = 1 3 = 1 3 = 1 

3 3 3 

E »'2;3-|-i-3;2,l V,„:3j E l'2:3 + i - 3;2 .2 V,„:3 j E ^"2:3 + 1 - 3;2 .3 V„;3j 

3 = 1 3 = 1 3 = 1 



E l'2;j + i-3;3,lV„-3j 
3 = 1 



E '"2;3 + i-3;3,2 V,„;3j E '"2 ;3 + i - 3 ; 3 ,3 V^, ; 3 j 
3 = 1 3 = 1 



3'"X3" 



Therefore, as in [5|, it can be proven that 



(2.10) 



^(m+l) X 



n{B)\ 



v: 



m+l I 



2.3. Periodic Patterns. This subsection studies periodic patterns in detail. 

For m,n > 1, a global pattern u — {aij)ij^z, on is called (m, n)-periodic if 
every i, j G Z, 



(2.11) 



for all p,q ^ Z. 

Let VB{'m,n) be the set of all (m, n)-periodic patterns and i?- admissible pat- 
terns. Let TB{m,n) — |PB(TO,n)| be the number of all (m, n)-periodic and B- 
admissible patterns. 

As in [3], VBijn^n) can be expressed by trace operators as follows. 

From (|2.9p . the periodic patterns in Y„i+i are given by Y,„_|_i:i, i ~ 1,5,9. 
Define 



(2.12) 



4=1,5,9 



m+l;i • 



Tm is called the trace operator of order m, as in [3]. Therefore, the following result 
is obtained. 

Proposition 2.1. Given B C S2x2(3), for m,n > 1, 
(2.13) TB{m,n)^tr{T'^J. 

Proof. The proof is similar to that for corner coloring in [2]. The details of the 
proof are omitted. □ 

Notably, from Proposition 12.11 V{B) ^ if and only if rB{m,n) > for some 
m,n > 1. 
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Recall some notation and terms from matrix theory. A matrix A is called nilpo- 
tent if A'^ = for some A; > 1. The property "nilpotent" can be used to specify 
whether _B is a cycle generator or non-cycle generator. 

Proposition 2.2. Given a basic set B C S2x2(3), 

(i) B is a cycle generator if and only iJTm is not nilpotent for some m > 1. 

(ii) S(-B) ^ % if and only ifVm is nilpotent for some ni > 1. 

Proof. From (|2.13p of ProDOsiton l2.11 B is easily seen to be a cycle generator if and 
only if tr(TJJj) > for some m,n > 1. Therefore, (i) follows immediately. 
Similarly, from (|2.10p . (ii) follows. 

□ 



The following proposition provides an efficient method to check the nilpotent for 
non-negative matrix and can be easily proven. The proof is omitted. 

Proposition 2.3. Suppose A is a non-negative matrix. Then, A is nilpotent if and 
only if A can be reduced to a zero matrix by repeating the following process: if the 
i-th row (column) of A is a zero row, then the i-th colume (row) of A is replaced 
with a zero column. 



3. Main Result 



3.1. Periodic Pairs. This section firstly classifies all local patterns in {0, 1, 2}^ 
into three groups. 



First, the local pattern a ~ 




ao,ai, 02,0:3) is assigned a number by 



3 

(3.1) (p{(ao, 01,02,03)) = 1 + y^aj3-^. 

3=0 

Then, all 81 local patterns are listed in the following three groups Go, Gi and G2. 
Go = {1,11,21,31,41,51,61,71,81} 

Gi = {2,3,4,7,10,12,14,17,19,20,24,27,28,32,33,34,38,40, 

42, 44, 48, 49, 50, 54, 55, 58, 62, 63, 65, 68, 70, 72, 75, 78, 79, 80} 

G2 = {5,6,8,9,13,15,16,18,22,23,25,26,29,30,35,36,37,39, 

43, 45, 46, 47, 52, 53, 56, 57, 59, 60, 64, 66, 67, 69, 73, 74, 76, 77} 



8 



HUNG-HSUN CHEN, WEN-GUEI HU, DE-JAN LAI, AND SONG-SUN LIN* 




Figure 3.1 

Clearly, every tile in Go can generate a (1, l)-periodic pattern. Furthermore, ele- 
ments in Gi, i — 1,2, form periodic pairs as in Fig. 3.1: two tiles that are connected 
by a line can generate a (2, 2)-periodic pattern. More precisely, the diagrams in 
Fig. 3.1 can be interpreted as follows. 

Proposition 3.1. (i) Each tile e in Go can generate a periodic pattern by 
repetition of itself: {e} is then a minimum cycle generator and is (1,1)- 
periodic. 
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(ii) For each tile e in Gi, there exist exactly three tiles ei, 62, 63 £ Gi such that 
{e,ei} can form a periodic cycle, which is a {2, 2) -periodic pattern, and 
{e,ei} is a minimum cycle generator for 1 < i < 3. 

(iii) For each tile e in G2, there exists exactly one tile e' G G2 such that {e, e'} 
is a minimum cycle generator and is (2, 2)-periodic. 

These minimum cycle generators are the simplest. 

Remark 3.2. From Fig. 3.1, the set Gi U G2 with 72 tiles can be decomposed into 
36 disjoint sets that each consists of two tiles such that each set is a minimal cycle 
generator. Therefore, the number of the tiles of maximal non- cycle generators in 
A/'(3) is equal to or less than 36. Indeed, elements in each pairs can be carefully 
picked up, and the maximum non-cycle generator with 36 elements thus obtained; 
see Table A. 3. Moreover, from Proposition \2.2l for each element B in these eight 
equivalence classes, = can be verified. 

3.2. Algorithms. Before the developed algorithms are presented , some notation 
must be introduced. 

Definition 3.3. (i) For a set A, let V{A) be the power set of A. 

(ii) For B CP(S2x2 (3)), let 

[B] = {[B] I B(EM}. 

(iii) For [B] € [P(E2x2(3))]; let (B) be a fixed chosen element of[B]. 

(iv) Let A/"* (3) be the set of all maximal non-cycle generators that can not gen- 
erate a global pattern. Indeed, 

J^*{3)^{N eJ\f{3) I S(iV) = 0}. 

(v) For B C £2x2(3), let C{B) be the set of all minimal cycle generators that 
are contained in B. 

(vi) For B C ^2x2(3), let M*{B) be the set of all maximal non-cycle generators 
that can not generate a global pattern and are contained in B . 

Now, the main idea of the algorithms is introduced, as follows. 
Let 

N^2^\ 

P(S2x2(3)) ^{B.j I < J < AT - 1}, where Bo = 0, 

Initial state for C(3) : C/(0) = {0}, 

Initial state for J\f*{i) : A/';(0) = {0}, 

Initial state for the set of aperiodic sets : Ui{Q) — {%}. 
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Main Algorithm 

repeat 

J = J + 1 

if V{B,) + 0, 

C/(j)=Cz(j-l)U{i?,} 
else 

if = 0, 

A/';(.7)=AA;(.?-i)u{i?,} 

else 

end 
end 
until j = TV 1 

After the algorithm has been executed, if Ui{N — 1) = {0}, then Wang's conjec- 
ture holds for p = 3. The methods to achieve the goal are introduced below. 

(I) reduce the number of cases that must be considered in the computation, 
(II) construct efficient initial states for C(3) and A/"* (3), 

(III) construct an efficient process for determining whether or not V(Bj) — 
and E(Bj) = 0. 

With respect to (I), the decomposition ^2x2(3) = Gq U Gi U 6*2 is used to 
reduce the number of cases that must be considered in the computation. Clearly, if 
-S C S2x2(3) contains a tile e S Gq, then i? is a cycle generator. Now, in studying 
Wang's conjecture, only cases S C Gi U G2 have to be considered. 

Given S = U A2 with Ax G P(Gi) and A2 G P(G2), if Ax or A2 is a cycle 
generator, then B immediately satisfies (II. 2p . By (|2.2p . the cases i? C Gi UG2 that 
have to be considered can be further reduced to the cases in X 01 1': 

1= {Ax\J {A2) I AxeVx and [A2] £ [P2]} 

(3.2) 

= {B, I 1<J<1X|} 

and 

(3.3) I' = {{Ax) U A2 I [Ax] G [Vx] and A2 G V2} , 
where 

(3.4) Vj = {A G P(Gj) I A ^ G for any G G C(Gj)} 

for i — 1,2. For brevity, the proof is omitted. From Table 3.1, N' = \I\ k, 
1.35075 X 10^2 and \T\ w 1.38458 x lO^^^ Therefore, X is the better choice for 
reducing B C Gi U G2. Notably, N' <C |P(Gi U Ga)! = 2^^ ^ 4.7223? x lO^i; the 
reduction is considerable. Table A.l presents the details. 

With respect to (II), let Ui{Q) = {0}. The initial data for C(3) are given by 
the set C/(0) of all minimal cycle generators that are the subsets of Go, Gi, or G2. 
Indeed, 

(3.5) CM = y C(G,). 
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On the other hand, the mitial data for TV* (3) are given by 

(3.6) A/';(0) = {TV G M*{Gi U G2) : |TV| = 36} . 

From Remark 3.2, Aff{0) equals the set of all maximal non-cycle generators in 
Gi U G2 with 36 tiles. C/(0) and TV/(0) can be easily found using a computer 
program. See Table A. 2 and A. 3. 

With respect to (III), the flowchart, which is based on (I) and (II), is as follows. 



f Cj(j) =Cf(j-l) 

{ Af;(j) = AT/u - ] 
( UiU) ^UiU - 1) 



Ci(j) = Cj(j-i) ^ 
A^-;(j) = A/-;o - 1) 




CiU) = c,u - 1) 

A'/O) = AT/U - 1) U [Sj] 



f cj(j) = ci(j-i: 
<^ Af;o) = Af;(j- 

[ Mf(j) = WjO - 1 





Initial data; C/ ( j — 1), J\ff(j — 1) and Ui(j - 1) 
Cosidcr oacli Bj G I, 1 < j < A' 




E(Sj) 7^ and V(Bj) = 
C7 0)=Cj{i-l) 

AA;o)=Ar;{j-i) 

UlU) ^UiU - 1) U [i3jl 



Figure 3.2 

Remark 3.4. Suppose that the computation based on the flowchart has been com- 
pleted. LetC ^Ci{N' -1), J\f* =TV;(TV'-1) andU ^Ui{N' - I) . 

(i) If the setlA — {0}, then Wang's conjecture holds for p — 3; otherwise, every 
element in hi is an aperiodic set. 

(ii) It is easy to see that an element in C may be not a minimal cycle generator. 
However, C(3) can be obtained from C = IJ [G] by the following process. 



cec 



If Ci, C2 G C with Gi C G2, then G2 must be removed from C. Indeed, 
^(3) = 1^ G C I C does not contain any G' ^ C except itself^ . 
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(iii) In a manner similar to that for (ii), let J\f = [J [N]. Now, 

A/"* (3) eJ\f I N is not a proper subset of N' for all N' e J\f except itself^ . 

Moreover, if U = {$}, 7V(3) = A/'*(3). 

3.3. Main result. The computer program of Fig. 3.2 is written, and the compu- 
tation is completed in finite time. Indeed, the cases that consume the most time 
are those in which the numbers of tiles in i? C Gi U 6*2 are 18 and 19. These cases 
can be computed completely within a week. The main result is as follows. 

Theorem 3.5. The set hi is equal to {0}, and Wang's conjecture holds for p — 3. 

Remark 3.6. (i) C(3) andj\f{3) can be obtained and their numbers are listed 
in Table A.4-- 

(ii) The computational results reveal that the maximum orders m of Tm[B) 
and V„i(i?) in applying Proposition \2.2\ are m ~ 35 and m =^ 13, respec- 
tively. More precisely, for B' = {2,5,13,36,53,60,62,64,77}, T^^iB') is 
not nilpotent but TkiB') is nilpotent for 1 < fc < 34. On the other hand, for 
B" ^ {2, 4, 5, 6, 9, 13, 14, 16, 18, 27, 32, 39, 60, 67, 78, 79}, V^iB") is nilpo- 
tent and Vfe(i3") are not nilpotent for any 1 < k < 12. The analytic proof 
that these numbers are maximal is not available. A prior estimate of the 
upper bound of m does not exist. 

For completeness. Tables A. 4 and A. 5 give the numbers of minimal cycle gener- 
ators and maximal non-cycle generators. 

Appendices 
A.l 

Table A.l presents the numbers of Vj and [Dj], j = 1, 2. Denote by 
gj{k) = the number of Vj with k tiles 

and 

gj{k) = the number of [Dj] with k tiles 
for 1 < j < 2 and 1 < fc < 18. 
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k 


ai{k) 


gi{k) 


.92 (fc) 


.92 (fc) 


1 


36 


1 


36 


1 


2 


576 


8 


612 


8 


3 


5304 


31 


6504 


34 


4 


31032 


146 


47988 


219 


5 


122184 


475 


256320 


971 


6 


342204 


1290 


998136 


3692 


7 


711288 


2581 


2812752 


10043 


8 


1129896 


4092 


5771988 


20554 


9 


1397892 


5005 


8886612 


31338 


10 


1361448 


4903 


10558368 


37319 


11 


1047816 


3763 


9807336 


34539 


12 


(i35r)80 


2321 


7125612 


23253 


13 


300888 


1106 


4007484 


14203 


14 


109080 


423 


1708632 


6162 


15 


29304 


118 


533664 


1945 


16 


5508 


28 


115164 


453 


17 


648 


4 


15336 


65 


18 


36 


1 


948 


8 



Table A.l 



A. 2 

Table A. 2 presents the equivalence classes of the minimal cycle generators in Gi 
and G2- 



k 


[C] e [C(Gi)] with k tiles 


2 


[{2,10}] 




[{2,40}] 


3 


[{2,12,19}] 




[{2,12.49}] 




[{2,42.79}] 



Table A.2 (a) 
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k 


[CI G \C(G2)] with k tiles 


2 


[{5, 3711 


3 


[|5, 45, 73)1 


4 


[{5, 13,30,4611 

L L " " " J J 




[{5,16.30,7311 

L L " ' " J J 




[{5,16,39,74)1 




[{5,9,46,64>1 




[{5, 13,35,64)1 


5 


[{5, 13,30,43,74)1 




[{5,15,30,52,73)1 

LL*'*' JJ 




[{5,9,39,46,76)1 

LL"'" ' JJ 




[{5,6,35,52,64)1 




[{5,15,35,46,66)1 


6 


[{5,13,35,43,57, 73)] 




[{5,6,43,53,66,73)1 




[{5, 13,36,53,66,73)1 




[{5,9,39,52,67, 74)1 

LL"'" 57 JJ 




[{5,9,39,43,74,76)1 




[{5,13,30,43,47, 64)1 




[{5,6, 13,43,47,66)1 




[{5,9,13,25,39,74)1 

LL"' ' w^^^j. JJ 




[{5,9, 13,30,47,64)1 

LL'*'' 'JJ 




[{5,6,16,47, 57,64)1 

LL ' " ' JJ 




[{5,9, 16,53,66,74)1 




[{5,9, 13,39,53,74)1 

LL'*''* JJ 




[{5,6,13,30,52,73)1 

LL * * ' ^ — J JJ 




[{5,9, 15,43,60,74)1 




[{5,13.35,45.66,74)1 


7 


[{5. 9. 13. 30. 32. 64, 74)1 




[{5,9, 13,47, 52,57,64)1 

LL"" " " 7" JJ 




[{5,6,16,36,53,66,73}] 




[{5,9,16,39,47,69,76}] 




[{5,6,13,35,43,66,73)1 




[{5,6,16,35,39,47, 76)1 




[{5,9,13,30,52,56,64)1 




[{5,6,16,35,36,57, 73)1 




[{5,9,15,22,46,56,66}] 




[{5,15,25,35,45,64,74)1 


8 


[{5,9,13,26,35,43,57,74)1 




[{5,9,13,35,39,52,74, 76)1 




[{5,9,13,45,47,52,56,64}] 



Table A.2 (b) 



A.3 

Table A.3 shows the equivalence classes of maximal non-cycle generators with 
36 tiles. 
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1. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 25, 

26, 27, 32, 33, 34, 35, 36, 42, 43, 44, 45, 52, 53, 54, 62, 63, 72}] 


2. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 25, 
26, 27, 32, 33, 34, 35, 36, 42, 43, 44, 45, 53, 54, 60, 62, 63, 72}] 


3. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 25, 
26, 27, 32, 33, 34, 35, 36, 42, 43, 44, 45, 54, 60, 62, 63, 69, 72}] 


4. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 25, 
26, 27, 32, 33, 34, 35, 36, 42, 44, 45, 54, 59, 60, 62, 63, 69, 72}] 


5. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 26, 
27, 32, 33, 35, 36, 42, 45, 57, 58, 59, 60, 62, 63, 68, 69, 72, 78}] 


6. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 26, 
27, 32, 33, 35, 36, 42, 57, 58, 59, 60, 62, 63, 68, 69, 72, 77, 78}] 


7. 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 26, 
27, 32, 33, 36, 42, 57, 58, 59, 60, 62, 63, 67, 68, 69, 72, 77, 78}] 


[{2, 3, 4, 5, 6, 7, 8, 9, 12, 13, 14, 15, 16, 17, 18, 22, 23, 24, 27, 
32, 33, 36, 42, 45, 57, 58, 59, 60, 62, 63, 66, 67, 68, 69, 72, 78}] 



Table A.3 



A. 4 

Table A. 4 shows the numbers of C(3) and JV{3). Firstly, denote by 

C3(fc) = {B G C(3) : \B\ = k}, 
Msik) = {N e : \N\ = k}, 

CsAk) = {[B] e [C{3)] : = fc for all B' e [B]} , 
AA3,e(fc) = {[Af] G [Af{3)] : \N'\=kiov all N' € [N]} . 

36 36 

Clearly, from Proposition 3.1, C(3) = IJ Csik) and Af{S) = \J Afaik). Only the 

fe=l k=l 

cases for C3(fc) ^ and JV3{k) ^ are listed. 



k 


|C3(fc)| 


|C3,e(fc)| 


10 


2880 


10 


9 


84600 


301 


8 


305388 


1094 


7 


2()4384 


952 


6 


105012 


406 


5 


21060 


102 


4 


3672 


29 


3 


528 


8 


2 


72 


3 


1 


9 


1 



Table A.4.(a) 
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k 






36 


1296 


8 


34 


720 


3 


32 


1152 


4 


31 


3168 


11 


30 


576 


2 


29 


288 


1 


28 


3168 


12 


27 


3456 


12 


26 


6048 


21 


25 


5760 


20 


24 


5184 


18 


23 


6624 


23 


22 


8640 


30 


21 


12672 


44 


20 


20160 


70 


19 


35280 


123 


18 


50256 


175 


17 


90000 


313 


16 


93024 


324 


15 


108720 


379 


14 


120384 


422 


13 


148536 


522 


12 


163512 


576 


11 


157536 


556 


10 


1SG48U 


()57 


9 


133200 


483 


8 


42624 


156 


7 


2160 


9 



Table A.4 (b) 



Acknowledgments 

The authors want to thank Prof. Wen- Wei Lin for suggesting the use of the 
concept of nilpotence to identify cycle and non-cycle generators. 

References 

[1] J. C. Ban and S. S. Lin, Patterns generation and transition matrices in multi- dimensional 
lattice models, Discrete Contin. Dyn. Syst., 13 (2005), no. 3, pp. 637-658. 

[2] J. C. Ban, W. G. Hu, S. S. Lin and Y. H. Lin, Zeta functions for two-dimensional shifts 
of finite type , Memo. Amer. Math. See, to appear.. 

[3] J. C. Ban, S. S. Lin and Y. H. Lin, Patterns generation and spatial entropy in two dimen- 
sional lattice models, Asian J. Math., 11 (2007), pp. 497-534. 

[4] R. Berger, The undecidability of the domino problem, Memoirs Amcr. Matli. See, 66 (1966). 

[5] K. CULIK II, An aperiodic set of 13 Wang tiles. Discrete Mathematics, 160 (1996), pp. 245- 
251. 

[6] C. Goodman-Strauss, Can't Decide? Undecide!, Notice of the American Mathematical 
Society, 57 (2010), pp. 343-356. 



DECIDABILITY OF PLANE EDGE COLORING WITH THREE COLORS 



17 



[7] B. Grunbaum and G. C. Shephard, Tilings and Patterns, New York: W. H. Freeman, 
(1986). 

[8] W. G. Hu AND S. S. Lin, Nonemptiness problems of plane square tiling with two colors, 

Proc. Amer. Math. Soc., 139 (2010), pp. 1045-1059. 
[9] J. Kari, a small aperiodic set of Wang tiles, Discrete Mathematics, 160 (1996), pp. 259-264. 
[10] A. Lagae AND P. Dutre, An alternative for Wang tiles: colored edges versus colores corners, 

ACM Trans. Graphics, 25 (2006), no. 4, pp. 1442-1459. 
[11] A. Lagae , J. Kari and P. Dutre, Aperiodic sets of square tiles with colored comers, Report 

CW 460, Department of Computer Science, K.U. Leuven, Leuven, Belgium. Aug 2006. 
[12] R. Penrose, Bull. Inst. Math. AppL, 10 (1974), 266. 

[13] R. M. Robinson, Undecidability and nonperiodicity for tilings of the plane, Inventiones Math- 

ematicae, 12 (1971), pp. 177-209. 
[14] H. Wang, Proving theorems by pattern recognition- II, Bell System Tech. Journal, 40 (1961), 

pp. 1-41. 

Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, 
Taiwan 

E-mail address: hhchen . amOOgSnctu . edu . tw 

Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, 
Taiwan 

E-mail address: wghuSmail . nctu . edu . tw 

Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, 
Taiwan 

E-mail address: werre216asfe87dirkagmail.com 

Department of Applied Mathematics, National Chiao Tung University, Hsinchu 300, 
Taiwan 

E-mail address: sslinSmath.nctu.edu.tw 



